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Abstract 

Insofar as the detection of anti-neutrinos from nuclear reactors is concerned, the 
SNO-I- detector - a 1 kilo-tonne liquid scintillator detector that inherits the experimen- 
tal infrastructure from the recently finished SNO experiment - is expected to perform 
just as well as the KamLAND experiment. The most important difference between 
these experiments is the distribution of nuclear reactors: whereas KamLAND has 9 
nuclear reactor sites within 300 km with a flux-averaged baseline of about 180 km, 
SNO-I- has only 1 within 300 km, with an average baseline of « 750 km. As a result, 
the reactor anti-neutrino flux at SNO-I- is only about 1/5 that at KamLAND, and the 
ability of SNO-f to constrain the solar neutrino oscillation parameter is diminished by 
a factor of about ■\/l/5 = 1/2.2 relative to KamLAND. In spite of this, SNO-K has 
comparable sensitivity to Am^2 ^ KamLAND because the rate of change of the spec- 
tral distortion as a function of this parameter is much greater than for KamLAND. In 
this report, this advantage is examined quantitatively using a geometric approximation 
that makes clear how the shape from SNO-I- has more statistical power than that from 
KamLAND. This result then is confirmed by determining the sensitivity to Amf2 using 
an ensemble experiment technique. 

1 Introduction 

The SNO-I- detector is a 1-kton liquid scintillator detector that can detect reactor 
anti-neutrinos via inverse beta decay. The performance of SNO-I- in the detection of 
these anti-neutrinos is basically identical to that of KamLAND: the target masses are 
almost identical, the detection efficiency and detector live time fraction are almost the 
same, and the background noise levels are similar. The main difference between the 
two experiments is the distribution of the sources of anti-neutrinos (Fig. [T|). Kamioka 
is surrounded by two reactor sites less than 100 km away and nine within 300 km, with 
a flux-weighted average of about 180 km; in contrast, SNO-I- has no reactor sites within 
a distance of 100 km, and just one within 300 km, with a corresponding flux-weighted 
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average distance of « 750 km. As a result of these differences, the reactor anti-neutrino 
flux at SNO+ is about 1/5 of that at KamLAND. 




Figure 1: The distribution of nuclear reactor sites in horizontal coordinates centered at 
SNO-|- (left) and at KamLAND (right). The center represents the detector position, while 
each point represents a nuclear reactor site, which may contain multiple reactors. The 
radius from the center represents the distance from the reactor, in logarithmic scale. The 
concentric circles show distances of 100 km, 300 km, and so on out to 13,000 km. Kam- 
LAND has two sites within 100 km and nine within 300 km, while SNO-I- has none within 
100 km and just one within 300 km. 

In spite of the lower flux at SNO-I-, its sensitivity to the neutrino oscillation param- 
eter Am^2 foi' ^ given exposure (say, in units of 10'^^ proton-years) is almost the same 
as for KamLAND. The statistical disadvantage of SNO-I- is compensated by the more 
distinct oscillatory signature of the anti-neutrino energy spectrum. In general, for a 
nearly background-free measurement above E^, — 3.4 MeV (this energy cut is applied 
to avoid geo-neutrinos, accidentals, and the ^■^C(a,n)^^0 background), the sensitivity 
to Amf 2 depends on the total number of detected events N and a factor K that de- 
pends on the rate of change of the shape of the distortion in the energy spectrum as a 
function of Amf2- The fractional uncertainty of Amf2 is proportional to the following: 

S (Amj2) K 

9"^ ^ -1= (1) 

Amf2 ViV ^ ' 

The statistical disadvantage of SNO-I- compared to KamLAND is expressed as 

\/ ^KamLAND / NO+ ~ ~ 2.2. The purposc of this document is to show that 

KKamLANo/ KsNO+ Las approximately the same value, so that the sensitivity of these 
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two experiments to Amjj is almost the same. 



2 Basic Information Used in This Study 

The latest result from KamLAND [1 is used as a baseline for the comparison with 
SNO+. Specifically, the exposure is chosen to be 2.44 x lO'^^ proton- years, and the 
best-fit solar neutrino oscillation parameter values was chosen as: 



Amfa = 7.55 x 10"^ eV^ (2) 
sin2 26i0 = 0.928 (3) 



The mixing angle expressed in terms of tan^ Oq is: 



tan^ 00 = 0.577 (4) 

These values are not exactly the same as given in [1 , but they are well within the 
uncertainties. 

The choice of uncertainty to use is not unique, since several different analyses are 
presented in [T]. Since the aim of this study is to compare the statistical power of 
SNO-I- with that of KamLAND, the errors for the analysis using only the KamLAND 
data are used (in preference to that from the global analysis of data from KamLAND 
and all solar neutrino experiments). Combining the statistical and systematic errors 
in quadrature, the la uncertainty on the oscillation parameters are as follows: 

Amj^ = (7.55t°JJ) X 10-^ eV^ (5) 



tan2 00 = O.577±EJ:J^° (6) 

The la range of sin^ 29q based on the above is [0.881, 0.973]. These ranges were used 
as a guide in determining the parameter range to examine in this study. 

The location and power of nuclear reactors was obtained using the information 
provided by the International Nuclear Safety Center (INSC) of the Argonne National 
Laboratory (ANL) [2]. The information was obtained in 2006; for details about the list 
used in this study, see Appendix|A] It is less accurate than information provided by the 
reactor operators, such as that used by the KamLAND collaboration. In particular, 
absolute normalization of the anti-neutrino flux is rather inaccurate because the rated 
power is often different from the power during operation, and because the INSC list 
has no information about the duty cycle. In this study, we used a constant factor of 
0.8 to represent the average duty cycle for all reactors {i.e. the power of all reactors 
was scaled down by this factor). 

The anti-neutrino detection rate depends not only on the flux, but on the detection 
efficiency. This was tuned in our model so that the number of anti-neutrino events from 
our model matches that from KamLAND. Although KamLAND 's detection efficiency 
is energy dependent, we ignored this dependence. As discussed below, a detection 
efficiency of 95% for E,y > 3.4 MeV was found to give a good match between the model 
and KamLAND. 
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The above inaccuracy of the nuclear reactor model also affects the predicted shape 
of the anti-neutrino spectrum after oscillations, but to a lesser extent that the overall 
normalization of the flux. Below, we shall show that the spectrum shape from our 
model is quite similar to that from KamLAND. Since the goal of this study is to 
compare the sensitivity of the two experiments to Am^j i the approximation should be 
adequate as long as the basic features in the spectrum shape are reproduced. 

3 Tuning the Nuclear Reactor Model to KamLAND 
Data 

According to the latest KamLAND result [T], the expected number of events with 
> 1.8 MeV passing all selection cuts before oscillations is 2179. In comparison, the 
reactor model used here predicts 3345 events, assuming 100% duty cycle and detection 
efficiency. Taking the duty cycle to be 80%, this becomes 2676 events. Taking the ratio 
of 2179 to 2676, one has a tuned detection efficiency for E^, > 1.8 MeV of 0.81. This is 
a reasonable value in view of KamLAND 's detection efficiency shown in the top frame 
of Fig. (2] 

In the study presented here, we shall focus on the essentially background-free region 
of Ei^ > 3.4 MeV, where KamLAND's detection efficiency is mostly well above 90%. 
In order to tune our model, we shall use as a standard the number of events after 
oscillations in the bin Ep S [3.6,4.0] MeV {E^ G [4.4,4.8]) MeV, which is about 160. 
Fig. [3] shows the model spectrum where the number of events after oscillations in 
the reference bin was tuned to be approximately the same as for KamLAND. The 
detection efficiency in the model after tuning was 95%. The agreement between the 
model and KamLAND spectra is reasonable. There is some disagreement around 
E^ = 3.4 MeV because the energy dependence of the detection efficiency has been 
ignored in the model. But this level of disagreement should not significantly affect 
conclusions regarding the relative sensitivity to Am^2 ™ KamLAND and SNO-I-. 

4 The Shape of the Reactor Anti-neutrino Spectrum 

The sensitivity of a given experiment to the neutrino oscillation parameters depends 
on how quickly the event rate and the spectrum shape change as a function of the 
parameter values. This is illustrated for KamLAND and SNO-I- in Figs. |4] and |5] 
There are a few notable features about the two figures. 

• The event statistics in KamLAND is about 5 times as large as that in SNO-I- . 

• In KamLAND, when Am^ is changed by Icr, the histogram contents in many bins 
changes by somewhat less than Icr. The same holds when sin^ 26 is changed by 
la. 

• In SNO-I-, when Am^ is changed by Icr, the relative change in shape of the 
spectrum in many bins is much larger than is the case for KamLAND. On the 
other hand, the change in the spectrum with a Icr change in sin^ 26 is much 
smaller than in KamLAND. 
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Figure 2: KamLAND's reactor anti-neutrino spectrum, adopted from Fig.l of ^Ij. The 
horizontal axis indicates prompt energy, which is E,y — 0.8 MeV. The vertical line at Ep = 
2.6 MeV {Eu = 3.4 MeV) divides the data into an almost background-free reactor anti- 
neutrino region above this energy and the region below which is contaminated by various 
sources of background noise. 



The above indicate that, because the event rate is much smaller in SNO+ than in 
KamLAND, SNO+ is much less sensitive to the mixing angle. However, because the 
rate of change of spectrum shape as a function of Am^ in SNO+ is much larger than 
in KamLAND, SNO-I-, for a given number of events, is more sensitive to Am^. Of 
course, this superiority in sensitivity is, to some extent, canceled out by the fact that 
the event rate in SNO-I- is only 1/5 that in KamLAND. 

The advantage that SNO+ has over KamLAND in its sensitivity to Am^ can be 
expressed quantitatively using a quantity K, which is related to the fractional error in 
Am^ as follows: 

— — ^ = " • —= , (7) 
Am^ -Jn 
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Reactor Anti-neutrino Spectrum at KamLAND 
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Figure 3: The model reactor anti-neutrino spectrum before and after neutrino oscillations 
The oscillated spectrum used the best-fit solar neutrino oscillation parameters in Eqs. [2](3] 



where N is the number of observed events, is a constant that depends on how quickly 
the shape of the anti-neutrino energy spectrum changes as a function of Am^, and a 
is a constant that is the same for any experiment (it depends on how the likelihood 
function is defined). The quantity K is given by the following expression: 



^-y'''^ HE.; Ami) 

In the above integral, the function f{E^; Am?) is the reactor antineutrino energy spec- 
trum for Am^ (the mixing angle is assumed to be fixed such that sin^ 29 = 0.928). It 
is normalized so that the integral is equal to unity: 



dE, f{E,;Am')^l (9) 

The integral for \/K^ has a form resembling x^, which is not an accident because 
it is extracted from a log- likelihood function for the fit parameter Am^. In applying 
this formula, one takes the red or blue function in Fig. [6] or [7] and the black function 
and calculate the difference squared. This is divided by the contents of the black 
function. The result is the integrand, which is integrated over the energy range of the 
fit. The more different the two functions, the larger the integral, which implies that K 
is smaller. Since K is proportional to the parameter uncertainty, one can conclude that 
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Figure 4: The change in spectrum as a function of neutrino oscillation parameter values in 
KamLAND. In both frames, the histogram in thick black is for (Am^, sin^ 20) = (7.55 x 
10"^ eV^, 0.928). In the left frame, the blue and red histograms show the change in the 
spectrum when Am^ is changed by la, while the mixing angle is kept fixed. In the right 
frame, the mixing angle is changed by la while Am^ is kept fixed. The error bar shows 
the statistical error in each bin for the thick black histogram. 

the greater the difference in shape for a unit change in parameter value, the smaller 
the parameter error. 

The value of K for KamLAND and SNO+ is shown in Table [l] The quantity K_ is 
the if -factor for the case where the function at the central value of Am^ is compared 
to that at —1(7, while is that for +la. The final row in the table shows the ratio of 
the factors, KamLAND to SNO+. The factor at KamLAND is about 2.6 to 2.8 times 
larger than at SNO-I-. This is a quantitative measure of how much more power SNO-I- 
has compared to KamLAND in extracting Am^ due to the rate of change of the shape 
as a function of Am^. Specifically, the ratio of the fractional error of Am^ between 
KamLAND and SNO+ is (from Eqn. [7| 

['^(^"^')/(^^')] KamLAND 
[5(Am2)/(Am2)]sNO+ 



KamLAND ^^q^ 



^SNO- 
(2.6 ~ 2.8; 
2.6 ~ 2.8 
2^2 



/SNO+ 


(11) 


D / ^SNO+ 


V %:amLAND 


TVS 


(12) 




(13) 
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[Reactor Anti-neutrino Spectrum KamLAND | [Reactor Anti-neutrino Spectrum KamLAND | 
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Figure 5: Same as Fig. |4j but for SNO+. 





K. 


K+ 


KamLAND 


10.3 


7.93 


SNO+ 


3.97 


2.85 


KamLAND/SNO+ 


2.59 


2.78 



Table 1: The K factors for KamLAND, SNO+, and the ratio of the factors, KamLAND to 
SNO+. K- is the factor for the comparing the function at the central value of Am? to that 
at —la, while is that for +la. The integral was evaluated over Eiy G [3.4, 10.0] MeV. 

= 1.2-1.3 (14) 

This is a remarkable result: although KamLAND has five times the statistics as SNO+, 
the latter's sensitivity to Am^ is 20—30% better because of the advantage in the shape 
of the spectrum, as indicated by the iiT-factors. The argument based on the X-factor 
is an approximation, but the conclusion is supported by a rigorous determination of 
the parameter sensitivity, which is discussed in the next section. 

5 The Sensitivity of the Experiments to Am^2 

The sensitivity of the experiments to Am^2 "^^s established rigorously by performing 
an ensemble experiment on randomly generated data. The random data were gener- 
ated using the central value of the neutrino oscillation parameters: (Am^2;Sin^ 29) — 
(7.55 X 10^^ eV^, 0.928). This was compared to the predicted spectrum at some other 
parameter pair using the following log-likelihood function: 
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Figure 6: The normalized reactor anti-neutrino energy spectrum for KamLAND for Am^ 
at the central value (black) and at itlcr. The mixing angle is such that sin^ 29 = 0.928. 



C=^[m-ni- log /ij + log r (n^ + 1)] (15) 

1=1 

The likelihood function was evaluated in the anti-neutrino energy range G [3.4, 10.0] MeV 
with 33 bins (bin width = 0.2 MeV). Energy below 3.4 MeV was excluded to avoid 
irreducible background noise, which significantly diminishes the importance of the con- 
tribution of this energy region to the extraction of the oscillation parameters. The 
quantity fit is the number of events expected for the current grid point in the oscil- 
lation parameter space, while is the randomly generated number of events for the 
central parameter pair. The gamma function is just a convenient way to evaluate nj. 
A total of 400 points in the parameter space was examined (Fig.jsj). At each grid point, 
data were generated randomly 1000 times and the log-likelihood function evaluated. 
The average and RMS of the evaluations was taken to represent the log-likelihood value 
and uncertainty at that point. 

The result for KamLAND is shown in Fig. [9] which is the marginalized log-likelihood 
(times 2) as a function of the oscillation parameters. The la errors from the ensemble 
experiment is obtained by determining the distance from the minimum to the point 
where 2 AC — 1; KamLAND 's ±1ct statistical errors are indicated by the vertical 
dashed blue lines. The errors agree fairly well. The la error for Amfa for KamLAND 
according to the ensemble experiment is: 
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Figure 7: Same as Fig. [6j but for SNO+. 
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Fig. 10 shows the marginahzed hkehhood as a function of the oscillation parameters 
for SNO+. As expected, SNO+ constrains the mixing angle much less because the 
statistics is much smaller, but the curvature of the likelihood as a function of Am^j is 
almost the same as for KamLAND. The la limit for SNO+ is: 



{7.55tof7) X 10-5 (eV)2 



SNO- 



(17) 



In Eqn. [14] the analysis based on the /^-factor predicted that SNO+ would have 
Am^ sensitivity that is 20~30% better than in KamLAND in spite of having only 
1/5 the statistics. Here, we find from Eqns. 16 and 17 based on a rigorous likelihood 



analysis that the relative sensitivity is about the same (lower error) or about 20% 
better (upper error). 



6 Conclusion 

Insofar as reactor anti-neutrino detection is concerned, the SNO-I- detector is expected 
to perform similarly to KamLAND. The main difference between the two is the distri- 
bution of nuclear reactors, which provide the detectors with anti-neutrinos. KamLAND 
has a greater number of reactors within 300 km compared to SNO+, resulting in a flux 
that is about 5 times greater. This gives KamLAND a natural statistical advantage 
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Figure 8: The grid in the neutrino oscihation parameter space over which the sensitivity 
of the experiments to these parameters was determined. The red star in the center of the 
plot indicates the position for the assumed true value of the parameters. The mixing angle 
grid was chosen to be equally spaced in tan^ 9. The grid in sin^ 26 is non-linear, and is 
indicated in red below the grid points. 

of \/5 ~ 2.2 in determining the neutrino oscillation parameters for a given exposure 
to reactor anti-neutrinos. SNO-I-, however, has an advantage based on the fact that 
relatively few reactors contribute to the bulk of detected events, and the few reactors 
that contribute to the bulk have baseline distances such that the oscillatory signature 
interfere constructively, resulting in the anti-neutrino spectrum shape varying rapidly 
as a function of Amfj- The advantage could be made quantitative in an approxi- 
mate manner using the so-called i^-factor technique. This factor quantifies the degree 
of similarity of the shape of a function at two different parameter values: the more 
similar they are {i.e. the slower the shape changes with parameter value), the larger 
K is. Thus an experiment with a small K is more sensitive to the parameter of in- 
terest compared to that with large K. It was found that K for SNO-f was 2.6^2.8 
times smaller than in KamLAND, which more than compensates for the statistical 
disadvantage factor at SNO-I- of 2.2. Combining the effect of statistics and the shape, 
SNO-I- is expected to outperform KamLAND in determining Am^ by about 20^30%. 
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A more careful likelihood-based analysis confirmed this result: the sensitivity of SNO-I- 
to Am^2 "^^s about the same or better by 20% compared to that of KamLAND. 

A The Position and Thermal Power of Nuclear Re- 
actors 

A list of nuclear reactor information was compiled based on information provided by 
the International Nuclear Safety Center (INSC) of the Argonne National Laboratory 
(ANL) [2]. As of 2006, there were 607 nuclear reactors in the list, although only a 478 
of them were reported as operational. The total rated thermal power of the reactors is 
1.3 TW, while the total of the 478 operational ones is 1.1 TW. A table summarizing 
this information is available at the following web site: 

http : //owl . phy . queensu . ca/~guillian/SN0+/reactor/2007-05-05/ 
reactor _list/Table .html 

The power flux from these nuclear reactors is shown in Fig. |11[ The power flux at a 
given location on Earth is given by: 



where Pi is the power of reactor number i, and Di is the distance between this reactor 
and the given point on Earth. 

More information is available at the following web site: 

http : //owl . phy . queensu . ca/~guillian/SN0+/reactor/2007-05-05/ 
index . html 

In particular, the list of reactors is available in text file format in the following sites: 

http : //owl . phy . queensu . ca/~guillian/SN0+/reactor/2007-05-05/ 
reactor _list/reactor _list_2006-ll-14.txt 

http : / /owl . phy . queensu . ca/^guilliaii/SN0+/reactor/2007-05-05/ 
reactor_list /react or _list_2006-ll-14_consolidated.txt 

The first list gives the thermal power, operational status, and latitude and longitude 
of each reactor, while the latter gives information about only operational reactors. 
Also note that the first one is a reactor-by-reactor list, while the latter is a site- by-site 
list, meaning that in sites with several reactors, the powers were combined. When 
combining the powers, only operational powers were considered. The total number of 
reactor sites is 208, whereas the total number of reactors is 607. 
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Figure 9: The marginalized log-likelihood value (times 2) as a function of tan^ 9i2 (top) 
and Ami2 (bottom). The red points are the average log-likelihood value over 1000 deter- 
minations, while the error bars show the RMS spread. In the bottom frame, the black 
curve through the points is the best-fit 3'''^ order polynomial. The vertical blue line is 
the location of the minimum, while the dashed vertical blue lines show KamLAND's la 
statistical error. The la error in this study is the distance from the minimum at which 
2AjC = 1. 
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Figure 10: Same as Fig. |9) but for SNO+. 
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